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' Abstract. It is shown that the radial part of the hydrogen Hamiltonian factorizes as the 

product of two not mutually adjoint first order differential operators plus a complex constant 
e. The 1-susy approach is used to construct non-hermitian operators with hydrogen spectra. 
^ , Other non-hermitian Hamiltonians are shown to admit an extra 'complex energy' at e. New 

self-adjoint hydrogen-like Hamiltonians are also derived by using a 2~susy transformation 
with complex conjugate pairs e, e. 
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O ; 1 Introduction 

The supersymmetric Quantum Mechanics (SUSY QM) has grown on the factoriza- 
C I tion and intertwining methods ^ apphed to transform the physical Hamiltonians. 

^ ' It yields new exactly solvable potentials which are either strictly isospectral to the 

CT' initial one because of broken SUSY, or almost isospectral due to unbroken SUSY (see 

the recent reviews The higher order SUSY QM amended the conviction that 
^ . the excited states can not be used to generate non-singular susy partners OEJElin]- 

Some applications deal with singular [7], soliton-type |H], periodic [Hj and other 
potentials ^Oj- Of special interest is the confluent algorithm |H1E] for which the 
second order procedure is applied to add a single level at an arbitrary point of the 
energy axis. However, almost all the works on the subject make use of transforma- 
tion functions with the real factorization constants e and with the factor-operators 
being always mutually adjoint. 

The case of complex e has not been studied to the desired extent. Exceptional 
cases seem jT21linilIlI , where the 2-susy treatment with e G C is formulated to ob- 
tain either hermitian or non-hermitian susy partners of a given initial Hamiltonian. 
Indeed, one of our purposes is to show that the Hamiltonians H can be factorized as 
the product of two not mutually adjoint first order operators A, B, plus a complex 
constant even for self-adjoint H. The method is not limited to the Hamiltonians 
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possessing the lower spectral bound; it can also yield non-hermit ian Hamiltoni- 
ans with complex potentials. Although the typical Hamiltonians of QM are hermi- 
tian, non-hermitian ones appear in Molecular Physics and Quantum Chemistry ^H] , 
Superconductivity [THj, Quantum Field Theory [T2j and other domains [T5]. The 
fact that they admit real eigenvalues for which the associated eigenfunctions are 
square-integrable JH] has been the basis of recent studies on PT-symmetry pmi21j . 
pseudo-hermiticity [2211231 and diverse physical models 

In this paper we illustrate these facts by constructing the hermitian and non- 
hermitian susy partners of the radial part (Hi) of the hydrogen Hamiltonian. Even 
though our non-hermitian operators are not PT-invariant, we shall see that a class 
of them has real eigenvalues identical to the hydrogen energies. The reality of the 
spectrum in this case is due to the breaking of supersymmetry. 

It will be shown that there is another class of non-hermitian operators having an 
extra square-integrable eigenfunction associated with e G C. In this case, the 'com- 
plex energy' e arises from the unbroken supersymmetry and, up to now, does not have 
a well established physical meaning (but see [25J). Unlike the phase-equivalent com- 
plex potentials [HI, the new 'bound state' associated with e is nodeless. Moreover, 
in counterdistinction with the formalism of PT-symmetry and pseudo-hermiticity, 
where complex energies appear in conjugate pairs, it turns out that e does not belong 
to the spectrum of the susy partner of Hi generated through e. 

In general, we shall see that the susy transformation is adequate to analytically 
determine normalizable eigenfunctions of non-hermitian Hamiltonians, including 
those with the complex energies. In this sense, our 'complex susy transformation' 
seems an analytical complement of the numerical techniques previously reported 
(see Pl] and references therein.) The eigenfunction connected with e is then removed 
by iterating the procedure in order to construct hermitian 2-susy partners of H^. 

The paper is organized as follows: Section 2 introduces the atypical factorizations 
Hi = AB + e, where e G C and the first order differential operators A and B are 
not mutually adjoint. Sections 3 and 4 are devoted to the construction of non- 
hermitian 1-susy and hermitian 2-susy partners, H((^) and H respectively, of H^. 
Final remarks and discussion are given in Section 5. 

2 The complex-type factorization method 

Let us consider a single electron in the field produced by a nucleus with Z protons. 
We shall use £ = Z/2rB and tb = /Ze^m for the units of energy and coordinates 
respectively. The corresponding time-independent Schrodinger equation reduces to 
Hii{)n,e{r) = EnilJn,i{r), with solutions 

V^n,£(r) =C„,,/+^ e-^/%Pi(^+ 1-71,2^ + 2; 2r/n), K = -1/n' (1) 
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where N 3 n = £ + s + 1; £ = 0, 1,2, ...,n — 1; s G Cn,i is the normahzation 
constant, i-Fi(a, c, 2) is the Kummer's function and //^(M"*", An) 3 ipn^iir) = rRn/{r), 
with an inner product defined by = 47r J^°° 'ijj{r)(f){r)dr < 00 and boundary 

conditions at r = 0: ip{0) = 0, tp'lr) = R{0). The effective potential Ve{r) has the 
domain Vy = [0, 00) and 

^.--^ + ^Kr) = -^ + V^--. (2) 

The Hamihonian Q is factorized as follows: 

He = AB + e (3) 

where the factorization constant is a complex number C 9 e := ei+ie2', ei, £2 7^ G M 
and the first order operators A,B, are not mutually adjoint (compare pi^): 

A:=-|- + /3(r), B:=^^+f3ir) (4) 

with f3 a complex-valued function fulfilling 

-/3'{r)+/3\r) + e = Viir). (5) 

The Riccati equation is solved by means of the logarithmic transformation (3{r) = 
— ^lnM(r) for which u is the most general eigenfunction of He (not necessarily 
normalizable) belonging to e = — /c^; C 3 k = ki + ik2; ki, k2 G M: 

u(r) = r^+^e~^^ f(r): 

JKJ, (6) 

/(r) := a iFi(£ + 1 - 1/k, 2i + 2, 2kr) + CU{i + l- 1/k, 2i + 2, 2kr) 

where a and ( are complex constants and U{a, c, z) is the logarithmic hypergeometric 
function. The global behaviour of these w-functions is analized in Appendix. 

Hence, for the /?-f unction we have: 

f^fr) = + k + nir); fi(r) := In /(r). (7) 
r dr 

A convenient expression for (3i{r), /32(r), the real and imaginary parts of (3{r) re- 
spectively, can be found in Appendix. 



3 New complex hydrogen-like potentials 

Let us consider the value e fixed. By convenience we shall make explicit the depen- 
dence of Q on (. Now, let us reverse the order of the factors in 0: 

BA + e = ~-^ + \4+i(r) + 2n'{r; C) = + V(r; C) := H(C) (8) 
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where H(^) is a non-hermitian second order differential operator and we have used 
(0) and ((Tj). The next step is to solve the related eigenvalue equation: 



H(C)^ = A^, A = Ai + iA2, Ai,A2 gM. (9) 

The dependence of \1/ and A on ^ will be dropped for simplicity. Notice that equa- 
tions (jni) and (jHl) imply an intertwining between the Hamiltonian Hi and the non- 
hermitian operator H((^): 

H(C)fi = M,, HiA = AV\{Q (10) 

Thereby, one sees that ^ oc Bif is a solution of Q if satisfies Hiip = Xip, while A 
reverses the action of B. Now, the general form of is obtained by taking e = — fc^ 
for A = — K^, and u{r) for (f{r) in (jHl). Hence, we have: 

ip{r) = r^+^e-^''{C iFi{i+l - 1/k, 2£ + 2, 2Kr) + DU{i+l - 1/k, 2i + 2, 2Kr)} (11) 

with C and D arbitrary complex constants. Therefore 

^ocB^ = ^^. (12) 
u 

We are looking now for the constraints on a, (, A, C and D leading to square- 
integrable First, consider A 7^ e {i.e., k, 7^ k); the behaviour of \1/ near the 
origin is 

I ^ r(m-i/K) ^ C 7^ 0, arbitrary a 

Thus, \l/ becomes divergent at r = except if either (I) A = — is real and 
K'^ = £ + s + l, sG Z^, or (II) A is complex but = 0. Let us pay some attention 
to these conditions 

CASE I (real A): Let us fix M 9 = £ + s + 1 = n G N. In this case, in fITTll 
iFi(a, c; z) and f/(a, c; 2;) are essentially the same function. Therefore, one can take 
D = and C = Cn/ (see equation ((l}), so ip{r) = ipn,e{r) and the functions (fT2|l 
behave asymptotically as \l/(r) oc ipn/ {f')r^oo- On the other hand, a straightforward 
calculation shows that these functions obey the following boundary conditions at 
r = 0: ^(0,0 = i?(0); ^'(0,C = 0) = -6eo; ^'(0,C ^ 0) = -£(^ + 5,i), with 
the kronecker delta. Thus 



^(r;C) oc 



1 , /iFi(£+l-n,2e + 2;2r/n) 

k h — In ' ^ ' ' 



n dr \ f(r) 



Mr) (14) 



are square-integrable eigenfunctions of H(C) with the real eigenvalues A = — 1/n^ 

En. 
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CASE II (complex A): For k G C and D = the function ()12|) behaves asymp- 
totically as follows: 



O (2^)^+1+1/-. r(^+l-l/«;) rV- ftl ^ U, <. ^ Q 



^ r(2^+2) e-'"- ^ , n / ^1 > 

(-2k)<+i-i/« r(£+l+l/K) r-1/- ^ U, <. ^ Q 



(15) 



which always diverges for r — >■ oo. Hence, there is no function * G L2(R+,47r) 
solving for a complex A 7^ e (z.e. k, ^ k) and Z) = 0. 

The formula (|T^ therefore gives all square integrable solutions {\['(r;(^)} of 
for A 7^ e. Concerning the case A = e, we see from equation (jH)) that any \E'e(r) in 
the one-dimensional kernel of A is an eigenfunction of H((^) belonging to e. After a 
simple calculation one gets 

^e(r) (X (16) 
u[r) 

This function can be in L^(M"'",47r) for appropriate values of a and (. In such a 
case, \E'e is an extra square- integrable eigenfunction of H((^) associated with e G 
C. Finally, a straighforward calculation shows that, although {\E'(r; (^), \E'e(r)} are 
elements of //^(M"*", 4-7?), they do not form an orthogonal set. A discussion on this 
kind of properties of the inner product is given in [21] . The next subsections analize 
these conditions and classify the resulting potentials according to their spectra. 



3.1 The real spectrum 

From Figure la, one sees that the behaviour of the complex function V(r;C), for 
a 7^ and C = 0, is given by 



V(r;C = 0) 




for r ~ 0, 
r ^ 00. 



Consistently with the case (A) of Appendix, the function \l'e(r) diverges at r = 0. 
Hence, \E'e(r) is not in L^(M+,47r) and the discrete spectrum ^"^(H^+i) of H(C = 0) 
is exactly the same as that of the hydrogen atom ad{\^e+i) = ad{Hi). Figure lb 
depicts the behaviour of tpn/ and the related function (jl4j) for one of the excited 
states. 

Another class of complex potentials V(r; () sharing the same spectrum as is 
obtained by considering the case (B) of Appendix. 
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(a) 




Figure 1: (a) The potential Ve+i{r) (dashed curve) and the real part of V(r;C) ior £ = 1, 

e = —(0.1 + iO.5)^, a = 1, and C = 0. (b) The corresponding unnornializcd |V-'3.i(^)|^ (dashed 
curve) and its susy partner |\E'(r;^ = 0)^ for E3 = —1/9 and the same values of the parameters. 

3.2 Complex potentials admitting complex 'energies' 

Let us consider ( ^ 0. Then the complex function V(r; Q behaves as shown on 
Figure 2a, i.e.: 

j Vf_i(r), for r ~ 0, arbitrary a 
V(r;C 7^ 0) ~ <( 

0, r — > 00. 





Figure 2: (a) The potential Vi_i(r) (dashed curve) and the real part of V(r;C), for i = 1, 



-(0.1 + z 0.5)2 



1, and C = 0.5 + iO.5. (b) The corresponding unnormalized |^'e(r)| for 



the same values of the parameters 

The case (D) of Appendix shows that, for (a, Q e Co and 9e{k) ^ 0, the function 
^e(r) behaves as 



( (2fc)^^+i r(^-i-i-i/fc) 
C r(2^-i-i) 



r ~ 



^g(r) < 



V — j_ — i^^^iK ^ KT r ^ 00, ki > O 

j.-1/k ^kr r ^ 00, ki < 



a r(2^-|-2) 
(2fc)^+l-Vfe 



(17) 



which clearly belongs to L^(M+,47r). Figure 2b depicts the global behaviour of 
|\E'ep. Therefore, the discrete spectrum (Td(H^_i) of H{( 7^ 0) is given by a^i^e-i) — 
ad{He) U {e}, provided that {a, () e Co, ee{k) ^ 0. 
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Finally, the case (C) of Appendix, gives another solution in the same class. 



4 The new real hydrogen-like potentials 

We shall now extend the analysis of the previous section by intertwining Hf^ with a 
new (to be determined) Hamiltonian H as follows 

HA = AHi (18) 

where the differential operator A is of second order 

and H reads 

H:=-^^+V{r). (20) 

The operators (|19H2m) depend implicitly on the label i. A straightforward calcula- 
tion allows to express the functions rj and 7 of ()19|) in terms of the auto-Backlund 
transformation of the solutions of (0) for ea and et |6j : 

Vir) = - f^TT^^^lTTl > + e,, (21) 



./3a(r)-A(r)^ 

7(r) = fi[{r) - Pl{r) + r^{r) (3,{r), ^ e, (22) 

Thus, the second order intertwining operator A in (fTU|) is expressed by two different 
solutions of the first order case. Moreover, it factorizes as A = a2ai, where 

a, = + = B; a2 = -^ + v-Pa = -A + r]. (23) 
ar ar 

Thereby, it is easy to rewrite A as 

A = {-A + ri)B = -Hi + e + riB (24) 

and potential is obtained through 

V{r) = Veir) + 2ri'{r) (25) 

In order to get a real ri{r) in (j2TJ|2SI), we consider the solution /3a(r) of (0) for G C 
as given and, by taking = and /?;,(r) = Pa{r) in (j2H), one finds 

Im(eJ _ €2 d 
Vir) = - ^ =-—-- = -—In uj{r) 26 

Im(/3„) f32{r) dr 
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where uj is defined in f)A.2|) of Appendix and tlie labels a and h have been dropped 
from e2 and /32- Henceforth, potential is real: 

V{r) = V,{r)-2^\nuj{r). (27) 

We are looking for potentials V{r) defined in the same initial domain Vy = [0, oo) 
(the situation when the initial domain is changed requires a different treatment, see 
e.g. Marquez et. al. in j7].) According to the proposition of Appendix, uj has at most 
one isolated zero in Vy. By choosing a proper u, the function u can be so constructed 
that its isolated zero coincides with one of the edges of Vy (compare [T^): 

limM(r) = =^ limco'(r) = (28) 

) — >0 r^O 



lim u{r) = ^ hm uj{r) = (29) 



r— >oo 



Let us examine the consequences. First, condition (j28|) is satisfied if the m in (jHl) 
is chosen with a ^ and ( = 0, so that the potential ^T7\ behaves as shown on 
Figure 3, i.e.: 

. ye+2ii") for r ~ 0, arbitrary fci 
l/(r;a^O;C = 0)~<! \. ' (30) 

r — > oo, arbitrary ki 



10 12 14 



(a) 




10 12 14 



(b) 



Figure 3: The initial potential Ve{r) (dashed curve) and its 2-susy partner V{r) for a = I, ( = 
and (a) ^ = 1, e = -(0.01 + i)^ (b) £ = 0, e = -(0.5 + i 0.1)^ 



On the other hand, condition ()29|) can be achieved for a = 0, ( ^ and ki > 0. In 
this case the equation (j77j) leads to 



V{r;a = 0-C 7^0) 



Ve{r) for r ~ 0, 
r ^ oo. 



(31) 



The real-valued potentials (jHOHjTl) resemble the hydrogen one and they could, in 
principle, represent physical systems as the Hamiltonians (|2Uj) are self-adjoint. The 
next step is to analize the new eigenvalue equation 



(32) 



8 



whose solutions, by similar arguments as for the first order case, are now obtained 
from the linear second order transformation (see equation ()18|)): 

^{r) oc AiJnAr) = {-En + e) + r]{r, Q ^(r, C) (33) 

where we have used fl24|) and \E'(r, Q is given by (fT^ . The corresponding boundary 
conditions at r = can be obtained from those of ipn/ir) and \E'(r, (). 

It is clear that the first term at the r.h.s. of ()33|) is in //^(R"*", An). The behaviour 
of the related second term is found by observing that 

1) If a 7^ and C = then rj diverges as at the origin while it is constant at 
r = oo. Therefore the product ri{r, 0) \E'(r, 0) is zero at the edges of and remains 
finite in all T>v- 

2) If a = and ( ^ then rj is a constant at both edges of T>v- Hence, 
rj{r, C 7^ 0) \I'(r, 7^ 0) is again well behaved in all T>v- 

Thus, the eigenfunctions ip oi H given in are in L^(R"'",47r) and H is an 
exactly solvable Hamiltonian with the same spectrum as the hydrogen atom. 

5 Summary and Discussion 

In this paper we have used a new type of factorization method to analyze a set of 
non-hermitian susy partners of the radial part of the hydrogen Hamiltonian. In order 
to generate the corresponding potentials we have used Darboux transformations 
with complex factorization constants. In contrast with PT-symmetry [T? H I ^ I^ 
and pseudo-supersymmetry [221, breaking of supersymmetry leads to purely 
real spectra. However, while the pseudo-hermiticity and PT-symmetry breaking 
involve pairs of conjugate complex eigenvalues, the unbroken supersymmetry (for 
which the non-hermitian 1-susy partners are not strictly isospectral) involves just 
a single 'complex energy'. To be more precise, in order to add two extra eigenvalues 
(real or complex) to a given spectrum one applies either twice the 1-susy procedure 
or a single 2-susy transformation (both can be made equivalent for the case we are 
dealing with.) Now, if the two new energies form a complex conjugate pair and if 
the functions ri{r) and 7(r) of the non-singular intertwining operator (fTI?|) are real, 
then the final Hamiltonian becomes self-adjoint and does not admit any complex 
eigenvalue, just as we have shown in Section 4. 

The problem of finding normalizable solutions belonging to complex eigenvalues 
for non-hermitian Schrodinger equations has been solved previously by numerical 
techniques in ^3] and analized inside a Lie-algebraic framework in [2ZI. In this 
paper we performed an analytical study and we expect that our results complement 
the numerical ones. Similarly as in other non-hermitian cases discussed in the 
literature, the interpretation of the 'complex energies' is an open problem, though 
notice possible applications to the absorptive (dissipative) systems j2Sl- 
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In general, we have shown that, the extension of the SUSY treatment to include 
complex factorization constants leads to results which are out of the scope of the 
PT-symmetry and pseudo-hermiticity. Indeed, the reality of the spectrum of the 
Hamiltonians H(C) in Section 3 depends on the parameter (, as it has been estab- 
lished in Subsections 3.1 and 3.2. As the non-hermitian Hamiltonians H{( ^ 0) 
in Subsection 3.2 do not satisfy the theorems by Mostafazadeh they are not 
pseudo-hermitian. On the other hand, it is not yet clear wheter the Hamiltonians 
H((^ = 0) of Subsection 3.1 could be pseudo-hermitian or not (though a primary im- 
pression can be depicted by noticing that the non-hermiticity of H(C) depends not 
on but on the non-trivial imaginary part of the factorization constant e.) Work 
in this direction is in progress and will be published elsewhere. 

Acknowledgements. The support of CONACyT (Mexico), project 40888-F, is 
acknowledged. 



Appendix 



The global behavour of the eigenfunctions uir) of equation © can be described in 
terms of a, C and the sign of fci: 

A) If a 7^ and C = then uir) is zero at the origin r = while it diverges at 
r = oo. 

B) If a = 0, C 7^ and k\ > 0, then u{r) diverges at the origin and tends to zero 
for r —>■ oo. 

C) If a = 0, C 7^ and ki < 0, then u{r) diverges at r = and r = oo. 

D) Let Co C C X C be the subset of complex pairs a 7^ 0, C 7^ 0, such that 

. , iFi(£ + l-l/fc,2£ + 2,2A:ro) w ^ ^ 
Uii + l-l/k,2i + 2,2kr,) ^ VroG(0,oo). 

^-(^)-^ + V([firiA) -""^ (A.1) 

is different from zero for (a, C) G Co , then u{r) is free of zeros in all T>v and diverges 
at r = and r = 00. 

E) If (a, C) G Co, ki < and Oi^k) = 0, then u{r) diverges at the origin while 
limr^oou{r) = (coinciding indeed with B). 

The presence of zeros in these functions has been studied by means of the following 

Proposition: Let u{r) G C^{D^) be solution of the Schrodinger equation 
u"{r) = [f (r) — e] u{r), where v{r) is a real- valued potential with domain Vy 



10 



and e G C. Assume that is a simply connected region of M. If Im (e) 7^ 0, 
then the complex- valued function u{r) admits at most one isolated zero in V^. 

Proof: Let 

:= (A.2) 
^ ' 2ilm(e) ^ ' 

where the bar denotes complex conjugation and W{-^ ■) corresponds to the 
Wronskian of the involved functions. Clearly uj is continuous on and uj'{r) = 
\u{r)\'^ > Vr G T>^, so uj{r) is always non-decreasing and can have either 
only one isolated zero or an entire interval of zero points in As every zero 
of u{r) is, necessarily, a zero of uj{r) then u{r) admits at most one isolated 
zero there, o 



The real function uj(r) in ()A.2|) plays a relevant role in the 2-susy approach of 
Section 4. A convenient expression for /5(r) in terms of u is given by 



Id,,,, f d ^ ^ 



/5(r)=A(r)+2/32(r) = --— lnu;'(r)+2e2(— lna;(r)) . (A.3) 
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